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Abstract 



We suggest a formula interpolating between the known asymptotic regimes of the BFKL 
equation as the approximate solution of that equation. The parameters appearing in this 
interpolation are fitted to the data on deep inelastic scattering in a wide range of the kinematical 
variables. Care is taken of the large-x domain as well, outside the HERA kinematical region. 
The boundaries and the interface between various dynamical regimes are also studied. 
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1 INTRODUCTION 

The Balitski-Fadin-Kuraev-Lipatov (BFKL) equation is a relativistic bound state equa- 
tion for two reggeized gluons. Although first applications concerned 77 scattering, the equation 
originally was derived for the hadronic (on mass shell) scattering amplitude in quantum chro- 
modynamics (QCD) at high energies yfs and fixed momentum transfer yf—i in the leading 
logarithmic approximation (LLA), implying the collection of all terms of the type (a s lns) n , a s 
being the QCD ("running") coupling constant. In this approximation, the total cross section 
increases rapidly, as 

°«Ut A=> I 1 - 1 ) 

V In s 

where u is (in the LLA) the position of the rightmost singularity in the complex angular mo- 
mentum plane of the t-channel partial wave with the vacuum quantum numbers (Pomeranchuk 
singularity), given by 

q 2 

uj = ^-Nln2 (1.2) 

TV 2 



for the gauge group SU(N) (N=3 in QCD), with a gauge coupling constant g = y/4ira s . 

Simply stated, the hope is that by solving the BFKL equation, one will be able to restrict 
the freedom available in the Regge pole theory, namely to determine the form of the leading 
singularities in the angular momentum plane, the form and the values of the parameters of 
the vacuum (Pomeranchuk) trajectory etc... Of all these ambitious expectations only (1.1) was 
of practical use, although still subject of uncertainties coming from the convergence condition 
a s In s < 1 and kinematical limitations. Any extrapolation of a s to t = should have discour- 
aged the direct application of (1.1) to fit measured hadronic cross sections. Really, even with 
a conservative value of u > 1.3, optimistic attempts in fitting data in the late 70-ies and early 
80-ies were soon abandoned. 

With the appearence of large virtualities, a new "hard" scale at HERA offered new possibil- 
ities for the interpretation of the "QCD-Pomeron" : most of the HERA results were claimed to 
confirm the existence of a " hard" Pomeron, manifest in the rapid rise of the structure functions, 
photoproduction of heavy vector mesons etc. . . While the experimental results were claimed to 
confirm quantitatively the predictions of the BFKL equation (rapid rise in s or in a; -1 = s/Q 2 ) 
and of the perturbative QCD in the whole, the interpretation proved once again to be mis- 
leading. Actually, the corrections to the LLA in the next-to-leading logaritmic approximaton 
(NLLA) were found large and distructive, lowering substantially the value of the Pomeron 
intercept |2j and raising doubts about the "hardness" of diffraction. 

What remains now from the " QCD Pomeron" is that it has a complicated j-plane structure 
and the rightmost (leading) singularity is located somewhere above unity. Since there is little 
hope that the whole perturbative series will ever be summed, and waiting for possible numerical 
solutions of the BFKL equation, all we can do now is to start from a "supercritical" (a(0) > 1) 
Pomeron and adjust it to the data. 

An important point to be mentioned here is that the above Pomeron is meant at the Born 
level, i.e. it should be subject to a subsequent unitarization procedure. Such a procedure is 
very complicated already in case of hadronic (on mass shell) reactions and becomes even more 
tricky as Q 2 7^ 0. A practical way out from this situation may be along the lines suggested in 
0, namely by introducing a Q 2 -dependent "effective" trajectory of a single, simple, factorizable 
supercritical Pomeron pole, absorbing the effect of its complicated j-plane structure as well as 
those from the possible unitarity effects (that cannot be calculated exactly anyway) in the form 
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of its "effective" trajectory. The free parameters are then adjusted to the experimental data. 
All what remains now from QCD is that a(0) > 1, which reflects the present status of the 
solutions of the BFKL equation at finite Q 2 . 

However, the solution of the BFKL equation is known exactly in the asymptotic, Q 2 — > oo 
limit as well, namely 

F 2 (x, Q 2 ) ~ exp y Jj 1 £n(l/x)eninQ 2 . (1.3) 

We suggest to use the two known solutions of the BFKL equation, (1.1) and (1.3) as boundary 
conditions and interpolate between these two by employing a minimal number of additional 
parameters. In our opinion this is the simplest solution of the problem, yet applicable to 
realistic processes at arbitrary values of Q 2 . By fitting the parameters to the experimental 
data, we hope to be consistent with unitarity and find the boundaries of different dynamical 
regimes, namely those governed by the BFKL equation or QCD evolution. 

In our previous paper Q we have already investigated such an interpolation within the 
GLAP evolution equation by assuming the "BFKL Pomeron" (1.1) to be the input, subject 
of a subsequent evolution. Actually, the asymptotic solutions of both the BFKL and GLAP 
equations have the same form (1.3) but the relevant paths may be different. Aiming at a high 
quality and reliable fit at small x, we have substanially improved as compare to our previous 
paper J|] the quality and range of the fits in the ("subsidiary") large-x domain. Still, the real 
paths of the solutions remain ambiguous for two main reasons: uncertainties in the low Q 2 , 
nonperturbative behavior (remaining outside of both the GLAP and BFKL equations) and 
unitarization (to effect both solutions), whose role increases with increasing Q 2 and decreasing 
x. Further work in this direction is needed. 

Our paper is organized as follows. In Sect. 2 we present our interpolating solution, in 
Sec. 3 we find the value of the adjustable parameters by fitting to the data on deep inelastic 
scattering; our conclusions are given in Sect. 4. 

2 THE MODEL 

We use the standard kinematic variables to describe deep inelastic scattering (DIS) : 

e(k) + p(P) -> e(k') + X , (2.1) 

where k,k',P are the four- momenta of the incident electron, scattered electron and incident 
proton. Q 2 is the negative squared four-momentum transfer carried by the virtual exchanged 
photon (virtuality) 

Q 2 = -q 2 = -{k-k'f . (2.2) 

x is the Bjorken variable 

x = 

W is the center of mass energy of the ("f*,p) 

W 2 = Q 2 ^—^ + m 2 (2.4) 

x F 

with m p being the proton mass. 

2.1 Structure function for low x and all Q 2 



-- . (2.3) 

2P.q K ' 

system, related to the above variables by 
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According to 0, we adopt the following ansatz for the small- x singlet part (labelled by the 
upper index S, 0) of the proton structure function, interpolating between the "soft" (1.1) and 
"hard" (1.3) regimes 

Ff 0) (x,Q 2 )= G(Q 2 )e A ^ 2 \ (2.5) 

with 



1+A(Q 2 ) 



and 



Ql 



(2.7) 



/_ i\/(<3 2 ) 
A(x,Q 2 )=[A(Q 2 )£n-j , (2.8) 

/(Q 2 ) = ^(l + e- Q2/Q ?). (2.9) 

The function f(Q 2 ) has been introduced in order to provide for the transition from the Regge 
behavior, where f{Q 2 ) = 1, to the asymptotic solution (1.3) of the BFKL evolution equation 
where f(Q 2 ) = 1/2. Alternative choices for this function, satisfying the boundary conditions, 
cannot be excluded, but we find our way of interpolation via (2.9) to be the simplest possible. 

It is customary to define an "effective" Pomeron intercept oup , which is in gereral x— and 
Q 2 — dependent, by rewriting the proton SF, introducing an "effective power" A e ^(x, Q 2 ) 

F^°\x,Q 2 ) = G 1 (Q 2 ) x- Ac// (^ 2 ) , (2.10) 

where the two effective quantities satisfy 

A e "(x,Q 2 ) =a e v ff (x,Q 2 )-l . (2.11) 

This definition with our parametrization leads to the identification 

G.iQ 2 ) = G{Q 2 ) , A^(x,Q 2 ) = ^fp- . (2.12) 

X 

It is worth noting that, in the limit 

/(Q 2 )~l, i.e. when Q 2 < Q\ , (2.13) 
the proton singlet (Pomeron component) SF reduces to 

Fi Sfi \x,Q 2 « Ql) ~ G(Q 2 ) x~^ . (2.14) 

We recover the standard (Pomeron-dominated) Regge behavior (with a Q 2 -dependence in the 
effective Pomeron intercept). Consequently, within this approximation 

A e "(x,Q 2 ) ~ A(Q 2 ) . (2.15) 

Therefore, at small and moderate values of Q 2 (to be specified from the fits (see below), the 
exponent A{Q 2 ) may be interpreted as a Q 2 — dependent (and of course a:— independent) effective 
power. 
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By construction, the model (singlet component) has the following Q 2 limits : 
- a) Q 2 — > oo, fixed x: 



F^°\x,Q 2 -> oo) -> A exp (JliMn^ £n^ ] , (2.16) 



which is the asymptotic solution of the BFKL and GLAP evolution equation (see Sect.l). 
- b) Q 2 -> 0: 

/n2 \ 1+A(Q 2 ^0) 

F f g2 ^ Q) ^ e A(,, Q ^o) V_ (217) 



with 



whence 



a 

A(g 2 ^0)^ e + 7 i72^ - e, (2.18) 

/(Q2 _ o) -> 1, (2.19) 

12 \ !+ e 



Ff°)(x,Q 2 -0)-A (£f (^j 6 ocQ 2 ^0 



^2.20) 



as required by gauge invariance. 

Apart from the (singlet, or Pomeron) component, a non-singlet component (sub- leading, 
or secondary Reggeons in terms of the Regge pole model) is also present at small x. Their 
contribution will be lumped in an "effective Reggeon" term, labelled by (NS,0) with the 
p, u>, f and A 2 reggeons absorbed in an effective trajectory with an intercept a r , namely: 

Fi NS >°\x,Q 2 ) = H(Q 2 ) X 1 "- , H(Q 2 ) = B ■ (2-21) 

The resulting low-x proton SF becomes 

Ff (*, Q 2 ) = F^°\x, Q 2 ) + Fi NS >°\x, Q 2 ) . (2.22) 

2.2 Extension to large x 

Our aim is to fix the free parameters appearing in our approximate solution of the BFKL 
equation (2.22). To this end, we could, in principle, fit (2.22) to the small-x data, where 
other contributions are negligeable. Whatever attractive, such a straighforward approach is 
not feasible since the small x domain where contributions to the SF other than (2.5) may be 
neglected is very narrow, making any fit unreliable. Moreover, the relevant limits depend on the 
models used, as shown explicitely in ||. Therefore we extend our model (2.22) by a "large- a;" 
part in order to have reasonable fits in a wide range of x. 

In our previous paper f|] the relatively simple and efficient model of Capella et al. [||] 
(CKMT) was used for that purpose. Below, aiming a better fit to the large- x data, we extend 
that model, by introducing additional adjustable parameters as follows 

F 2 (x,Q 2 )= F^°\x,Q 2 )-(l-x) p ^+ F 2 {NS >°\x,Q 2 )-(l-x) R ^ . (2.23) 
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Similar to we use the following Q 2 -dependent exponents of the large x factors (with 6 
additional parameters) 



P(Q 2 



Poo + 



PQ-Poo 

i + QVQl 



R{Q 2 



Too + 



1 + Q 2 /Q 2 



(2.24) 



2.3 Total cross-section for (7*,p) scattering 

The structure function is related to the total cross-section of virtual Compton scattering 
(or approximatively to the transverse cross-section, if the longitudinal component is neglected) 
by 

4n 2 a 1 /. 4m 2 x 2 ' 



'tot 1 



x,Q* 



1 + 



F 2 (x,Q' 



<f(W,Q 2 ) 



(2.25) 



Q 2 1 - x V Q 2 , 

where x is replaced by using (2.4). 

In the limit Q 2 —>■ 0, the SF vanishes as Q 2 and we write the total cross-section for (7,p) 
scattering (with real photons), as a function of the center of mass energy W 



rto?(W) = 4n 2 a lim Q2 



F 2 (x = Q 2 /W' 2 ,Q 
Q 2 



A-n 2 a (A a' 1 '' W' 2e + B b~ ar W' 2 ^^ 



[2.26) 



with W 



12 



w 2 



ni 



3 FITTING TO THE DATA; RESULTS 

In our fitting procedure the experimental data sets from "HI" 



"E665" 0, "NMC" H, SLACJTI 



B0,"ZEUS'' [0,0,0, 

'BCDMS" [0 for the proton structure function F 2 {x, Q 2 ) 
were used as well as data points [M on the (7,p) total cross-section <Jt2f\w), in the kinematical 
region with < Q 2 < 5000 GeV 2 , < x < 0.75, W > 3 GeV. The total number of data for 
this "complete" ensemble is 1253 (see Table 1). 
The resulting fits may be commented as follows. 

The use in || of an "economic" (only 8 free parameters) large-x extension a la CKMTQ] 
resulted in good fits {Xdof ~ 0.67) within a rather limited domain in x < 0.1, with a selection 
of about 310 data points (mainly those of the low x HI data). 

Better fits at large x (including the BCDMS and SLAC data) can be achieved (at the 
price of 6 additional parameters with the large x extension, presented above. Namely, we 
obtained Xd.o.f ~ 1-16, distributed among each subset of data as shown in Table 1. This is a 
considerable improvement with respect to the small x, all Q 2 results of [|]]. Interestingly, the 
present improvement of the large-x behavior has little effect on the values of the parameters 
governing the low x behavior and fitted previously |4]]. As noted, 6 new parameters were 
introduced in our present version of the model (see Table 2). Among a total of 16 parameters, 
14 are free (instead of 8 in H), the remaining 2 being fixed in the following way: 

1. Similar to [|], we choose the "canonical" value |1| e = 0.08; 

2. Similar to [|J, we estimate from QCD the parameter 71 = 16iV c /(ll — 2//3) with four 
flavours (/ = 4) and three colors (iV c = 3), it equals 5.76. It corresponds to the asymptotic 
regime (when Q 2 — > 00, or f(Q 2 ) — > 1/2)), far away from the region of the fits, where / = 1, 
hence the value 71 = V5.76 = 2.4 is more appropriate in the domain under consideration. Fits 
with 71 free confirm this choice. 
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To compare with, the proton structure function and (j*,p) cross section are calculated in 
the ALLM model in the whole experimentally investigated kinematical range. The ALLM 
fits to the data are good (Xd.o.f ~ 1-1, recalculated) with a total of 23 adjustable parameters. 
This is also to be compared with the 21 parameters used in 0, where a quite good fit (with 
only a slightly better Xd.o.f) * s a l so presented. 

The results of our fits for the structure function versus x for fixed selected Q 2 are shown in 
Figs. 1-2. 

The total cross section for real photons on protons as function of W 2 is displayed in Fig. 3. 

The agrement with the data is impressive except may be for the SF reported at the high- 
est Q 2 - value (=5000 GeV 2 ). We recall once again that this result is obtained with 14 free 
parameters. 

Q— slope as a function of x 

Known as a useful tool in the sudies of low x region, the derivative of the SF with respect 
to InQ 2 

d / n2 N dF 2 {x,Q 2 ) 

(Q— slope for brevity) measures the amount of the scaling violation and eventually shows the 



transition from soft to hard dynamics. It was recently extracted from the HERA data [13]. 
where the variables x and Q 2 are strongly correlated, because it is implied that, for a limited 
acceptance (as it is the case in the HERA experiments) and for a fixed energy, one always has 
a limited band in Q 2 at any given x, with average < Q 2 > becoming smaller for smaller x. 

From the theoretical point of view however, this derivative depends on two variables (x 
and Q 2 ) which are quite independent and one is not restricted to follow a particular path 
on the surface representing the Q— slope. That is why we discussed in || the Q-slope as a 
three-dimensional quantity. Here we adopt a more pragmatic attitude and we only compare 
in Fig. 4 the predictions of our model with the available experimental results for the special 
experimental set of (x, Q 2 ) chosen in | I3| . Our predictions are quite in agreement with the 



data, as it may be expected from a good fit of the structure function implying also agreement 
with its experimental derivatives : an increasing Q— slope when x decreases down to ~ 10~ 4 , 
then a "turn over" : the slope decreasing at lower x (or Q 2 ) values. Notice that the turn over 
point, located at Q 2 ~ 2 GeV 2 , may be related to the different behaviors of the gluon and qq 
sea distributions. 

x— slope as a function of Q 2 

The derivative of the logarithm of the SF with respect to Enl/x 

d(m(l/x)) 

(a;— slope for brevity), when measured in the Regge region, can be related (for low x) to the 
Pomeron intercept. In Fig. 5, we compare the calculated a;— slope with the quantity called 
"effective power" (A e // ). In fact, the measured (x— dependent) quantity should not be confused 
with the "effective power" introduced in (2.10-12) (A e ^(x, Q 2 ) in our notation), it is really the 
x— slope. The agreement between experiment and theory is very good. 

Finally, in Fig. 6 the Q 2 -dependence of this calculated derivative together with the power 
A (2.8) is shown for some low x - values. On the same figure, the behavior of the function 
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f(Q 2 ) (2.9) is also shown. In our model, Regge pole behavior is equivalent to the condition 
that f{Q 2 ) is close to unity 

This lower limit, marked on Fig. 6 (tentatively approximated within a 2 % accuracy for the 
function f(Q 2 )), is located near 60 GeV 2 . Until this landmark, the power A indeed remains 
very close to the x— slope. Beyond, Regge pole behavior is not valid (since / ^ 1) and A 
cannot be considered as the x— slope any more. On the other hand, the x— slope turns down 
as Q 2 increases, approaching its "initial value" of ~ 0.1 at largest Q 2 and coming closer to the 
unitarity bound. Notably, at large Q 2 the derivative gets smaller as x decreases, contrary to the 
general belief that dynamics becomes harder for smaller x, but in accord with an observation 



made in pl| . Care should be however taken in interpreting the "hardness" of the effective 
power outside the Regge region. 

The transition region occurs when f(Q 2 ) goes from 1 to 1/2 i.e. in a band in Q 2 estimated 
between ~ 60 GeV 2 and ~ 5000 GeV 2 . 

4 CONCLUSIONS 

The consequences of the present studies are manifold. 

On the one hand, the present approximate solution of the BFKL equation ((1.5-9) with 
the parameters listed in Table 2) can be used in the future to compare it with the numerical 
solutions of this equation as soon as they will be available (and there is a strong need for such 
numerical solutions !). 

Further, our interpolating formula may clarify the still open question: what is the Pomeron ? 
While the answer is known within the context of the analytic S'-matrix theory (that gave rise 
to the notion of the "Pomeron"), namely that it is an isolated moving j-plane singularity with 
vacuum quantum numbers and with a(0) = 1, an alternative definition may arise from QCD, 
namely that the Pomeron is the solution of the BFKL equation. In that case the deviation from 
the simple s a<yt '® > behavior may be indicative of the difference between the two. The model 
presented in this paper shows it. 

Another aspect of the Pomeron studies is its "hardness" at low (or vanishing) Q 2 . It is 
well known (see e.g. [22]) that in the present energy range, the cross sections (or structure 



functions) can well be described by logarithmic functions as well (rather than powers in s (or 
1/x)). Apart from numerical fits, this phenomenon has also a physical interpretation P^ ]: in 
the presently available energy range there is sufficient phase space available only for a finite 
(small) number of gluons rungs in the Pomeron ladder, each contributing with a power in 
£ns (or £n(l/x)), resulting thus in a £n(s), £n 2 (s), ... behavior of the cross sections (same 
for the SF). As repeatedly stressed, such logarithmic parametrizations, being equally efficient 
phenomenologically, have the advantage of being consistent with the unitarity bounds. Note 
that the extrapolation between the low Q 2 logarithmic behavior and high Q 2 asymptotics may 
be more complicated than that presented above. For that purpose the model of pi| , combining 
the logarithmic behavior at Q 2 = with a "supercritical" power behavior far off shall may be 
appropriate. 

Finally, we remind that the " BFKL-Pomeron" presented in the present paper is not a direct 
solution of the BFKL equation. Instead, it interpolates betwen the two known asymptotic 
solutions of that equation and as such is a small step forwards with respect to those previously 
known. Our poor knowledge of other solutions (in different kinematical regions ?) of the BFKL 
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equation shows on one hand the complexity of the problem, but on the other hand it is clear 
that further progress in theory cannot be achieved without a better understanding of such 
basic objects in high energy physics as the bound state of two gluons or the elastic scattering 
amplitude. 

Acknowlegdments. We thank V. Fadin, E. Kuraev and L. Lipatov for numerous discus- 
sions on the Pomeron. 
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Observable 






Number of points 


x 2 


x 2 


Experiment — year of pub., Ref 






i 

L|_JI 


Ft 












HI -95 [| 


1 




93 (0 in m 


77 


_ 


#1-96 [R 


1 




193 


149 


110 


HI -97 [10 




44 


55 


20 


ZEUS - 96 


11 




188 (0 in |§) 


256 


— 


ZEUS - 97 


12 


i 


34 (0 in U) 


13 


— 


ZEUS - 98 


13 




44 (0 in H) 


28 


— 


£665 - 96 \] 


L4| 




91 (0 in i) 


101 




NMC -95 [ 




156 (0 in §) 


177 




SLAC - 90/92 




136 (0 in |U) 


145 




BCDMS - 89 




175 (0 in g) 


257 




*lo! [18| 


99 (73 in g) 


177 


73 


Total 






1253 (310 in §) 


1434 


203 


X 2 I d.o.f. 








1.16 


0.67 



Table 1. 

Observables used in the fitting procedure. The complete available experimental kinematical 
range is taken into account for Q 2 {i.e. 0. < Q 2 ( GeV 2 ) < 5000, x (2. 10~ 6 < x < 0.75), with 
the lower limit W > 3 GeV. 

Also shown is the distribution of the partial \ 2 f° r each subset of data used in our fit with 
the parameters listed in Table 2. The results from §1 are also shown. 



low X 

A = .1365 
a (GeV) 2 = .2372 

7 2 = .02158 
Q2 (GeV) 2 = .1991 
Qf (GeV) 2 = 1490. 

B = 1.944 
b (GeV) 2 = 1.804 
a r = .3603 
large x 
Q2 (GeV) 2 = .1864 
Po = -26.32 
Poo = 10.52 
Q 2 (GeV) 2 = 13.86 
r = 2.601 
= 3.846 



Table 2. 

Free parameters used in our fit. The non-fitted values are e = 0.08 (fixed from JOJ) and 
7i = 2.4 (suggested by QCD). 
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Fig. 1 Proton structure function F 2 (x,Q 2 ) as a function of x at low Q 2 . The data shown 
are listed in Table 1, the error bars represent the statistical and systematic errors added in 
quadrature, the solid curves are the results of our fit (the parameters are listed in Table 2). 
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Fig. 2 Proton structure function ^(x, Q 2 ) as a function of x at various values of Q 2 (see also 



Fig. 1). The dashed lines are the results of the ALLM model [20 . 
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Fig. 3. Limiting case to the real photon-proton total cross-section o^t as a function of W 2 , 
square center of mass energy (see also Fig. 1). The dashed lines are calculations from the ALLM 
model j20|. 
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Fig. 4. Q-slope = dF Q^cp " as a function of x (for the indicated Q 2 values). The experimental 
points are from |13 |. The predictions in the present model are given in continuous line. 
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Fig. 5. A e //(< x >,Q 2 ) experimental points from as a function of Q 2 (for the indicated 

< x > values), compared with the calculated x— slope values = d Q^n(i/x)) • ^he predictions in 
the present model are given in the approximation including only the low x Pomeron contribution 
(2.5). 
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Fig. 6. x— slope = ^"in/vaO) as a f unc ^ on °f Q 2 (f° r the indicated x values). The predictions in 
the present model (solid line) are calculated in the low x Pomeron approximation (see Fig. 5). 
On the same left scale the exponent A (2.8), equivalent to the Pomeron intercept-1 when 
f(Q 2 ) ~ 1 is also plotted. This function f(Q 2 ) is also shown in dash line (right scale). The 
estimated transition region is between the vertical landmarks (see the text). 
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